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Abstract 

The emitted power of the radiation from a charged particle moving uniformly on a circle 
inside a cylindrical waveguide is considered. The expressions for the energy flux of the radia- 
tion passing through the waveguide cross-section are derived for both TE and TM waves. The 
results of the numerical evaluation are presented for the number of emitted quanta depending 
' on the waveguide radius, the radius of the charge rotation orbit and dielectric permittivity 

. of the filling medium. These results are compared with the corresponding quantities for the 

^'j^ ' synchrotron radiation in a homogeneous medium. 

^ ; 1 Introduction 

■ , The wide applications of the synchrotron radiation (see, for example ||l], |2[ ^] and the references 

Q I therein) motivate the importance of investigations for various mechanisms of control of the 

' radiation parameters. From this point of view, it is of interest to investigate the influence of 

a medium on the spectral and angular distributions of the synchrotron emission. This study 
I is also important with respect to some astrophysical problems H]. The radiation from a 

■ charged particle circulating in a homogeneous medium was considered by Tsytovich in Q (see 

^ ' also 0, I). It had been shown that the interference between the synchrotron and Cerenkov 

radiations leads to interesting effects. The radiation from a charge rotating around a dielectric 
ball enclosed by a homogeneous medium is investigated in papers |l^. The presence of the 
ball leads to interesting effects: if for substance of the ball and particle velocity the Cerenkov 
condition is satisfied, there are strong narrow peaks in the radiation intensity. At these peaks 
the radiated energy exceeds the corresponding quantity in the case of a homogeneous medium by 
some orders of magnitude. A similar problem for the case of the cylindrical symmetry we have 
considered in Refs. [0, |T^. In the paper a recurrent scheme is developed for constructing 
the Green function of the electromagnetic field for a medium consisting of an arbitrary number of 
coaxial cylindrical layers. The investigation of the radiation from a charged particle circulating 
around a dielectric cylinder immersed in a homogeneous medium [jl2[, has shown that under 
the Cerenkov condition for the material of the cylinder and the velocity of a charge there are 
narrow peaks in the angular distribution of the number of quanta emitted into the exterior space. 
For some values of parameters the density of the number of quanta in these peaks exceeds the 
corresponding quantity for the radiation in the vacuum by several orders. 

The present paper is devoted to the radiation from a charge rotating inside a cylindrical 
waveguide filled by a homogenous dielectric. The contribution of waves of the electric and 
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magnetic type to the total radiation intensity is determined. The results of the corresponding 
numerical calculations are presented and a comparison with the intensity of the synchrotron 
radiation in a homogenous medium is carried out. 



2 Electromagnetic fields inside a waveguide 

To determine the electromagnetic field of a charged particle rotating inside a cylindrical waveg- 
uide with conducting walls and filled by a material with dielectric permittivity £q , firstly we shall 
consider a more general problem. Let a charge q moves with the velocity v along a circle of radius 
po in the plane z = inside a cylinder of radius pi and axis z. The permittivity of the medium 
inside the cylinder is Eq and it is immersed in a homogeneous medium with permittivity £i ( for 
simplicity we shall assume that the permeability is equal to 1). The 4-potential Ai = {cp, —A) 
of the corresponding electromagnetic field can be determined via the Green function (being a 
second rank tensor) by the formula 

Mx) = -^^J Gii{x,x')ji{x')d'^x', x = {t,r), i,l = 0,l,2,3, (1) 

where a summation over index / is assumed, c is the speed of light in the vacuum, and ji{x) 
is the 4-vector of the current density of the charge. In the cylindrical coordinates {p, (/?, z) the 
spatial components of the latter are in form 

31 = — 5{p - pq)5{(p - ujQt)5{z)5i2, v = ljqpo, Z = 1,2,3 (2) 

PO 

(the values of indices / = 1,2,3 correspond to the coordinates p,yp,z respectively). A recurrent 
scheme for the construction of the Green function in a medium consisting of an arbitrary number 
of coaxial cylindrical layers is developed in [11|. The Green function for the problem under 
consideration can be derived from the general formulae of this paper. By using the problem 
symmetry we can write the following Fourier expansion for the Green function: 

+ 00 „ 

Gii{x,x')= ^ dkzdujGii{LO,m,kz,p,p')exp {ilmiip - If') + kz{z - z') - uj{t -t')]} . (3) 

m=— oo 

In the Lorentz gauge for the components we need here and for p = pQ < pi, p < pi from the 
general formulae of Ref. one has 

TT 



Gii{uJ,m, kz, p, Po) = ^7[Zi X! Jm+a{>^OP<)Hm+a{>^OP>) + B^m^Jm+a{>^OP) 

a=— 1,1 



/ = 1,2, 



(4) 

where /9< = min(p, po); P> = max(p, po)) Jm{x) is the Bessel function, Hm{x) = Hm\x) is the 
Hankel function of the first kind. In (Q) the coefficients Bm^ are defined by the expressions 

Tjia) T 

- -«^m+a(,AoPoj-f7777 f? \ + 

^ iaXlHm{XiPi)Hm+a{XiPi) Jm+p{XoPo) 

7rpl(3iW{Jm+a,Hm+a) ^^-^W {Jm+p, Hm+pY 

where, as in the paper we use notations 



c^kl 



-^JmiXoPl) K,,,, ^ rr (7) 



£1 — £0 2 ^^^^ W{Jm+l,Hm+l) 
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with 



W{a,b) = a{Xopi) — ^ o(Aipi)- 



(8) 



dpi dpi 

In formula @ the summand with m = does not depend on time and hence does not give a 
contribution to the radiation field. Therefore in the consideration of the radiation field we can 
assume m ^ 0. 

Substituting Eq. (|2[) and expressions (^), (^) for the Green function into formula (|l|) we 
receive 



-l-CXJ 

Mx)= J2 e^'"('^--°*W dk,e'''^'Ami{k,,p), 

m=-oo 

where in the Lorentz gauge 



(9) 



Aml{kz-,p) 
A.m?,{kz,p) 



vq 



T { ^Sl + E O^'B^^Um+aM \, 1 = 1,2, 
a=±l 



(10) 



Here the term 



A 



(0) 
ml 



a' Jm+a ( AoP< ) Hm+a ( AoP> ) , 



-1.1 



corresponds to the field of a charge in a homogeneous medium with permittivity Eq. 

The vector potential of the electromagnetic field inside a waveguide with conducting walls 
can be obtained from the formulae given above taking the limit ei — > cxd. From formulae @-(|8|) 
it follows that Ai(ei) enters only in the argument of the Hankel function. By a simple calculation 
in view of the asymptotic formulae of this function for large values of the argument [13| it can 
be seen that in the limit ei — > cxd the coefficients (^) take the form 



-Jm+a{XoPo) 
iCk-Jrn.— 



H. 



m+a 



(Ao/5i) 



Jm+a{XoPl) 



7rpiAoJm(AoPl)Jm(Ao/Ol) 



p=±l 



Jrn+p{XoPo) 
Jrn+p{XoPl) 



(12) 



Thus, the field inside a cylindrical waveguide is determined by formulae (^), (|10|), where the 
coefficients Bm^ are defined according to Eq. (12). Note, that in the limit po —>■ pi one has 
Bm^ = —Hm+a{Xopo), and in virtue of Eq. (pT|) fields (|lO|) tend to zero. We could expect 
this result as when the charge approaches the wall of a waveguide the particle and its mirror 
image cancel each other. Before to proceed further calculations we shall consider the analytical 
properties of Ami{kz,p) in Eq. (^), as a function on the complex variable k^- From formulae 
(^), (|TT|), (|l2|) it can seem, that the points kz = ^muo^/eo/c are branchpoints for this function. 
However, by using the expansions of the cylindrical functions (see, for example, [^]), it can 
be seen that actually A^iikzjp) is a function on Ag only: Ami = Amii^Q, p). Note that the 



logarithmic singularity In Aq contained in the expansion of the Hankel function in is cancelled 



with the corresponding singularity of the first term in Eq. (12). As a result the function 
Ami{kz,p) is meromorphic on the complex plane kz- From the formulae (|ll]), ( |T^ ) it follows 
that this function has poles corresponding to the zeros of the Bessel function Jm(AoPi) and its 
derivative: 



kz — i^l ^ 



2, ,2 



£0 



Jm.n 
2 ' 
Pi 



a 



n 



1,2, 



(13) 
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with the notation Jm\x) = d'^ Jm/dx^ . Note that the singularities in the first and second 
terms on the right of formula (p^) corresponding to the zeros of the functions Jm±i{^oPi) are 
cancehed, and consequently the functions Bm^^ are regular at these points. In Eq. ([T^ jm,n are 
positive zeros of the Bessel function (a = 0) and its derivative (cr = 1), arranged in ascending 
order, < Jmn+i- these zeros are simple, and hence the values of given by Eq. 

(13) correspond to simple poles of the function Ami{kz,p)- They describe the eigenmodes of a 
cylindrical waveguide and are known as TM modes in the case a = and TE modes in the case 
o- = 1 0. 



For real values of eo poles (13) are situated on the real axis of the complex plane k^, for 
jm,n < mwo/Oi-y/eo/c, and are purely imaginary otherwise. In formula (P) to derive unambiguous 
result for the integral over k^ it is necessary to specify the rules to escape the real poles (see 
Section |^. For this purpose note that in physically realistic situations the permittivity is 
complex: Eq = e'q + iEq, where the imaginary part Eq > describes the absorbtion in the 
medium. From here it follows that the radicals ^ki^^ from Eq. (13) with a positive/negative 
real part are situated in the upper /lower half of the complex plane k^. This leads to the following 
rule to circle the poles in (^): in the integral over kz the positive poles should be circled from 
below, and negative ones - from above. Note also that if the walls of the waveguide have 
finite conductivity, then to the wave number k^ an additional imaginary part is added (see, for 
example, |14]), which leads to the same rule to escape the poles. 



3 Radiation intensity 

In this section we shall consider the radiation field travelling inside the cylinder at large dis- 
tances from the charge. First of all let us show that in Eq. ([l^) the fields and the term 



corresponding to the first summand in definition (|12| ) do not contribute to the radiation field. 
This immediately follows from the estimate of the integral over kz in Eq. (^) by using the sta- 
tionary phase method (see, for example, |jl^). As in this integral the phase kzZ has no stationary 
points, for large \z\ the integral tends to zero faster than any power of 1/ \z\, provided that the 
coefficient of the exponential function belongs to the class C°°(R). It follows from here that the 
the singularities of this coefficient can only give the contribution to the radiation field. As it 
has been mentioned in the previous section, for the integral over kz in Eq. (^) the singularities 
correspond to the poles of the second summand in Eq. (12) at k^ = ifcn''', defined by relations 



(|T3|). To determine the corresponding contributions to the integral we note that, as it has been 
mentioned above, by taking into account the absorption, the points kz = kn\—kn'') are situ- 
ated in the upper (lower) half of the complex plane kz- Therefore the integration contour over 
kz in Eq. ^ can be closed by a semicircle of a large radius in the upper (lower) half-plane. This 
choice is stipulated by the fact that for large \z\ the integrand exponentially tends to zero in the 
upper (lower) half-plane for z > 0(< 0). As a result the corresponding integral vanishes when 
the radii of the semicircles go to infinity. Thus, for large \z\, z > (<)0 the integral over kz in (|9|), 
according to the residue theorem, is equal to the sum of the residues at poles kz = kn\—kn'') 
multiplied by 27ri. For large n, when jm\ > mwoPi-y/eo/c, the poles ^k^^ are purely imaginary 
in the limit eg ^ and the corresponding contribution tends to zero exponentially for z oo. 
As a result these poles do not contribute to the radiation field. Therefore, the radiation field far 
from the charge takes the form 

, (<^) 
Ai{T,t)=2^i e^-(^-°*) ^ E I^es^.=fcMA„^e*^^^ (14) 

m=—oo (T=0,1 n=l 
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where the maximal value nmax is defined by the conditions 



J M < V^om , 



.{a) V pi 

J (.) > V^m . 

ni,nmax+l C po 



(15) 



Thus, at large distances from the charge the radiation field inside the waveguide is presented 
in the form of waves with discrete set of values of the wave vector projection on the waveguide 
axis, 



1 2 r?^*^^ 

i , , . . . , ( (-max ) 



, defined by formula (jl^). Having Ami we can find the scalar 
potential by using the Lorentz gauge condition. As a result the electric and magnetic fields can 
be presented in the form of sums of the TE and TM waves: 



+00 



m=— oo 



(7=0,1 n=l 

For the z-components of the fields corresponding to the elementary waves one has 

2g Jmijm,nPo/pi) 



(16) 



E 



(0) 
m3 



Jm {£]nP/pi ) , H';^i = 0, for TM waves, 



r(0) 



H. 



(1) 

m3 



'Jm\jm,n 

'iqvi jm,n Jm {jm,np/pi ) {jm,nPo/pi 



E. 



0, for TE waves. 



^Pl kn'^ {Tn? jrn]n) '^mO'"^!") 

The transverse components can be found from the formulae (see, for example, []l^ ) 



(17) 
(18) 



E 



(0) 
mt 



H 



(1) 

mt 



ikn ^ Pi 
jm,n 

ikn ^ Pi 
jm,n 



V,*, H 



(0) 
mt 



V.*, E, 



,(1) 

mt 



ck^^ 



ck 



(1) 



^Sti-mt 



(19) 



where = E^^ for TM waves and ^ = H^^l for TE waves, Vt = {d/dp, im/ p, 0), and es is the 
unit vector along the axis z. The energy flux through the cross section of the waveguide per 
unit time interval is given by the Poynting's vector S: 

pi r^TT 



I 







pdp I d(p{e3 ■ S) 







S = -[EH1, 



(20) 



where the fields are defined by expansions ([iq). Substituting these expansions into Eq. (|20| ) and 
using the formulae for the integrals involving products of the Bessel functions [13] it can be seen 
that the contribution of the terms E^''(j^,L,p)H^,'' (j^/^^„/,p) is proportional to (5m,m"5rm'5o-cr'- 
As a result the intensity / can be presented as a sum of the radiation intensities on separate 
modes: 



JO) 



^ / J / J mn ' / J mn ' 



(21) 



m=l yn=l 



n=l 



where the energy radiated on the frequency w = muj^ in the form of the TE and TM modes per 
per unit time interval is determined by the formulae 



.(TM) 
/(TE) 



''max 



Er(TM) r(TM) 
'■mn 1 'mn 



2q^Vm kn^ JmUm,nPo/pi) 
"f-" Jm,n "^m+lw"*!™/ 

2g2f^m 



(0) 



n=l 

'^max 



V- r(TE) r(TE) 
/ y mn ' mn 



.(1)2 
Jm,n 



Jm{jm,nPo/Pl] 



"^PlPO k^rt\j^S - m2) Jlij^r^l) 



n=l 



(22) 



(23) 
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In these expressions the terms with fixed n correspond to the waves with a given value = 
of the wave vector projection on the waveguide axis. Note that in the hmit po — > pi the 



radiation goes to zero. Expanding the Bessel functions in the numerators of Eqs. (22) and ( P3|) 
in this hmit we receive 

4™> - (24) 

eo \ PlJ 

(TE) _ 2g^u;gm (jffn - m^) po^ 

^mn (I -) (25) 

under the assumption jm}n{l — po/ pi) <^ I- As a necessary condition for the presence of a 
radiation of a given type on the frequency to = mujQ one has 

^ (cr) (a) Jm,n . ^ 

nwo > uj'^ i, uj'^ n = c j=, (26) 

Pi VEO 



where we have introduced the boundary frequency oim.ri for given m, n and for a given type of 
waves {a = 0,1). In terms of the charge velocity it looks like u-^/eo/c > {jm\/m){po/ pi). In 

particular, by taking into account the inequality > m (see, for example, [jl3[), as a necessary 
condition for the presence of a radiation flux through the cross-section of the waveguide we have 
Vy/£o/c > po/ Pi- If this condition does not satisfied we have an interesting situation when the 

rotating particle does not radiate. Note that the boundary frequency uJm]n is a characteristic 
of the waveguide and does not depend on the parameters of the charge (energy, radius of the 
orbit). In formulae ( |2^ ) the quantities kn^ are expressed through the boundary frequency by 
the formula 



A;W = v^v/m2u;2-u;^l'. (27) 
From Eqs. (22), ([27|) it follows that if the charge is orbiting with frequency ujq = ujm]n/'rn, i.e. 



if the frequency of radiation coincides with the one of boundary frequencies (26), the radiation 



intensity for the TE waves goes to infinity. Note that in this limit the projection kz = kn^ goes 
to zero. However, it is necessary to take into account that under these conditions the absorption 
in the medium and in the walls of the waveguide becomes important, and consequently it is 
necessary to take into account the imaginary part Eq of the permittivity of the filling medium 
and finite conductivity of the waveguide walls (the damping constant in the waveguide walls 



depending on the conductivity is presented, for example, in |14| ). Thus, formulae (22) are valid 
for frequencies, not too close to the boundary ones, when kn^ » I3\, raiOQ^^l^j where /?a is the 
damping constant due to the ohmic losses in the walls of the waveguide. 

Introducing 7m,n = jrn^n)po/^Pi, the arguments of the Bessel functions in the numerators 



of Eqs. (22), ( ^3D can be written as m^m^n- Note that for these quantities one has the following 
inequalities Pq/ pi < ^m^n < v^/e^jc. Let us consider the radiation intensities (p2|), ( [2^ ) for large 



values of m. From the well known properties of the Bessel functions for large values of the order 
it follows that the behavior of these intensities essentially depends on the sign of 1 — ")m,n- 

First of all assume that the Cerenkov condition is not satisfied, Vy/eQ/c < 1, and synchrotron 
radiation is present only. For l — v^/£Q/ c ^ 1 the radiation intensity has maximum for c<j ~ o^c = 

mcWo, with nic = — '^ml%max)~'^^'^ i ^-nd as it follows from Eq. (|l5| ) nic < {l — v'^so/c^)~^^'^- For 
the frequencies lu ^ Uc the emitted power exponentially decreases. If the Cerenkov condition 
is satisfied, for the modes with 7m,n > 1 in addition to the synchrotron radiation we also have 
Cerenkov radiation. Note that in the waveguide the sufficient condition for this radiation at a 
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given harmonic m is Jm^nmax > 1 instead of v^/Io/c > 1. For this modes the radiation intensity 
for large m approximately linearly increases in the mean with frequency. This is a characteristic 
feature of the Cerenkov radiation. In this case to obtain a finite result for the total intensity 
( ^Tl) it is necessary to take into account the dispersion for the permittivity. Consider now the 
limiting case of large values of the waveguide radius /Oi — > oo. In this limit the main contribution 
to the radiation intensity is due to large values jm,n, when we can use the asymptotic formula 



m — 1 (—1) , 
n + — — + ^-^ . (28) 



From this relation it follows that AAq = (jmii+i — im,n)/pi = t^/pi- Now dividing and multi- 
plying expressions (22) by vr/pi and taking the limit pi ^ oo the sum over n can be replaced 



by the integration over Aq. As a result one has 



,(TM) q^vm p^oy^/c ^Jm^uj^eo/c^ ~ ^2 m (oo\ 
lom = / T J^(AoPo)dAo, (29) 

£0P0 Jo AO 



^Om = L / ^ . , ^ =^»/m(AoPo)^Ao. (30) 



c'po Jo 



^m^u^eo/c^ - A§ 



By introducing a new integration variable 9 according to Aq = {muJQ^/£Q/c) sin 6 for the total 
radiation intensity (including the regions z > and z < 0) with the frequency to = ttiujq in the 
limit pi — > oo we receive 

JOm — ^{^Om "T -'Om ) — „ /— ^ 

cv^o 



X 



7r/2 







cot^ 9 j!^{mv ^/eq sin 9 / c) H — ^ ^o<^m("^^\/^si'^^/'^) 



sin 6ld6' (31) 



This formula coincides with the expression for the radiation intensity of a point charge circulating 
in a homogeneous medium ^. 

We carried out numerical calculations for the number of the emitted quanta per one period 
of the particle orbiting, 

Y(TE) _ 2^ r(TE) jy(TM) _ 27r (TM) (o^) 

nmujQ hmujQ 

for various values of the parameters Eq, pi/po- In figures Q and || the quantities Nm^\ -/V?t7™^ 
for the harmonic m = 24 are presented as functions on the ratio pi/po for an electron with the 
energy 2MeV and for the media with permittivities Eq = 3 and Eq = 1, respectively. The narrow 
peaks in the graphs of the number of quanta for TE waves correspond to the singularities on 
the boundary frequencies Wm.n. For the values of parameters corresponding to figure || one has 
v^/eqIc ~ 1.67 and the Cerenkov condition is satisfied. As j^i/m ~ 1.24, j^i/fn ~ 1.1, then 



the condition ( |26D is satisfied for all values of the ratio pi/ po > 1. As a result the radiation is 
nonzero for any pi/po > 1- For values of the parameters corresponding to fig. 2 the Cerenkov's 
condition is not satisfied {v^/eq/c ~ 0.97) and the radiation corresponds to the synchrotron 
emission. Note that in this case the radiation is much more weaker in comparison to the case 
of figure ||. Now there is a range of values of the ratio pi/po for which condition ( p6| ) is not 
satisfied and in this region the radiation is absent. It is well seen in fig. |2|. We also carried out 
numerical calculations for an energy of an election equal to 0.6 MeV and for the same values of 
the remaining parameters, as in figure 1. In this case the Cerenkov condition is not satisfied and. 
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Figure 1: The number of quanta emitted on the harmonic m = 24 in the form of the TE and TM 
waves per circulating period of an electron, multiplied by fic/q^ , hcNm^'^ /q"^ and hcNm^^ /q^, 
versus the ratio pi/pQ- The energy of an electron is 2 MeV, and the dielectric permittivity for 
the filling medium is eo = 3. 




8 



as calculations have shown, the qualitative behaviour of the Nm and Nm in dependance of 
the ratio pi/po is similar to the case presented in figure |2|, and the radiation intensity is much 
more weaker to compared with the case of figure Jl]. 

For comparison we have presented in figure y the dependances on the permittivity of the 
number of quanta emitted from an electron in the form of the TM and TE waves in a waveguide 
and in a homogeneous medium on the harmonic m = 24. The energy of an electron is 2 MeV, 
and pi/po = Note that the intensities of the both types of radiation increase when the 




1.5 2 2.5 3 3.5 4 4.5 5 1.5 2 2.5 3 3.5 4 4.5 5 



Figure 3: The number of emitted quanta ficNm/q^ for the TM and TE waves on the harmonic 
m = 24 as a function on the dielectric permittivity eo in the waveguide and in the homogeneous 
medium (dashed lines). The energy of an electron is equal to 2 MeV, and pi/ po = 1.5. 

permittivity increases. 
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